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EPHEMERIS OF A HIGHLY ECCENTRIC ORBIT: EXPLORER 28 


Barbara E. Lowrey 


ABSTRACT 

An ephemeris has been obtained for Explorer 28 (IMP 3) which agrees well 
with 2 years of radio observations and with SAO observations a year later. This 
ephemeris is generated over the 3 year lifetime by a numerical integration 
method utilizing a set of initial conditions at launch and without requiring further 
differential correlation. Because highly eccentric orbits are difficult to com- 
pute with acceptable accuracy and because a long continuous arc has been ob- 
tained which compares with actual data to a known precision, this ephemeris 
may be used as a standard for computing highly eccentric orbits in the Earth- 
Moon system. 

Orbit improvement was used to obtain the initial conditions which generated 
the ephemeris. This improvement wa.-s based on correcting the energy by adjust- 
ing the semimajor axis to match computed times of perigee passage with the 
observed. This procedure may generate errors in semimajor axis to compen- 
sate for model errors in the energy; however this compensation error is also 
implicit in orbit determination itself. 
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EPHEMERIS OF A HIGHLY ECCENTRIC ORBIT: EXPLORER 28 


I. INTRODUCTION 

Highly eccentric orbits offer unusual opportunities for fitting long arcs with 
high precision. One major advantage is that a high eccentricity orbit is seldom, 
or never, affected by the earth's atmosphere. Since the atmospheric drag force 
is poorly predicted, the uncertainties in close earth circular orbits have to be 
assumed to arise from this source, and indeed, close earth orbits are often used 
to derive the density of the atmosphere. In the case of high eccentricity orbits, 
however, precise determinations can be obtained and extrapolated forward with 
good accuracy. 

A second major advantage is provided by the good resolution of the perigee 
position of the orbit. Because the raw data will yield a good definition of the 
position of perigee, the time of perigee passage becomes an important param- 
eter in the orbit improvement. This parameter has the additional benefit of 
being readily derived from both the reduced observational data and the numerical 
computation, thus facilitating comparisons between the two. Further, the in- 
vestigation of the long term behavior of this parameter yields some insights 
which may not be clearly visible in short arc fitting. 

IMP 3 (Explorer 28) satellite is a good example of a high eccentricity orbit. 
Radio tracking has been obtained for over 120 orbits, that is, 2 years of the 3 
year lifetime. This paper discusses a technique of improving the orbit by 
finding a long arc which is consistent with the available short arc data. An 
ephemeris has been obtained which gives a good fit to the available data (time 
of perigee passage). This ephemeris, when extrapolated ahead one year, proved 
capable of making predictions of accuracy such that SAO's Baker- Nunn network 
was able to acquire the satellite. The technique of obtaining a good fitting 
ephemeris is discussed in this paper. 

Earth satellites in highly eccentric orbits (e > .9) are subject to large and 
erratic perturbations by the moon, as well as a strong solar perturbation. Also, 
the earth's oblateness term may cause a serious computational error, even 
though the secular effect of the oblateness is quite small. It can be difficult to 
distinguish between a legitimate perturbation and a numerical error introduced 
by improper selection of integration interval, or a programming error. 

The IMP 3 (Explorer 28) orbit provides a good sample of the major per- 
turbational effects in the earth-moon system. Since there is a total of 2 years 
data with which a numerical integration has been compared, it is suggested that 
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the ephemeris of IMP 3 serve as a computational standard. For most applications 
the full accuracy of the ephemeris will not be required; however, a program com- 
parison with the IMP 3 ephemeris will test the ability of the program under test 
to compute an orbit with high eccentricity and large semimajor axis in the Earth- 
Moon system within desired accuracy requirements. 

The general history of the satellite is discussed in Section n. Section HI 
gives a summary of the perturbations, the characteristic changes of specific 
orbital elements and the origins of the changes. Section IV describes the method 
of computation and discusses the short period effect of the earth’s oblateness, 
which is the major source of computational errors and consumption of machine 
time. Section V explains the data which is employed in the comparison and in- 
dicates the technique involved in obtaining initial conditions which will give a 
satisfactory fit over the entire arc. Section VI gives the comparison between 
the data and the numerical integration. Section VII presents the ephemeris. 

The final section, Section VIE, presents the conclusions. 

II. HISTORY OF SATELLITE 

The satellite was launched on May 29, 1965. The initial orbital elements, as 
assumed in this report, are 

a Q = 21.726138 earth radii 

e 0 = 0.95251083 

i Q = 33?828446 

w Q - 135?72696 

fi 0 = 138?45267 

M 0 = 0?023663 

or, in vector form, 

x Q = 6099.5844 km 
y 0 = 602.05128 
z Q = 2409.1608 
x Q - 1.1047527 km/sec 
y 0 = 9.8556127 
z Q = -4.4520836 
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Epoch time is May 29, 1965 at 12:07 hours U.T. (a = semimajor axis, e = ec- 
centricity, i - inclination, w = argument of perigee, ft = right ascension of the 
ascending node, M = mean anomaly. The angular elements are taken with re- 
spect to the earth's equator.) The initial value of perigee is about 200 km above 
the earth’s surface, and the apogee is about 240,000 km. The orbital period is 
about 140 hours, or 5.8 days. 

Tracking data was acquired for the first two years of the lifetime. After 
the cession of radio transmission in March 1967, there were no observations 
of any kind until March 1968 when the Smithsonian Astrophysical Observatory 
was able to photograph the satellite with Baker-Nunn cameras (Figure 1). At 
this point, the perigee had been lowered enough to permit visual observations. 
These observations, which were taken during several perigee passages, on 
March 21, 27, April 8, May 7 and May 30, confirm the validity of the ephemeris 
during the unobserved period in the last year of the orbit. 


III. PERTURBATIONS 

The major perturbations on a highly eccentric satellite orbit are due to the 
gravitational fields of the sun and the moon. In addition the oblateness and higher 
order terms in the earth's gravitational field and the solar radiation pressure 
may produce a discernible effect on the ephemeris. 

The height of perigee is the element which is critical to the success of the 
mission. The lunar-solar perturbations on perigee height are large and must be 
carefully considered when selecting the launch conditions to ensure an adequate 
lifetime for the mission (Reference 1-4). Both the gravitational fields of the sun 
and the moon cause the long term trend, often termed the free oscillation, in 
the eccentricity and inclination as a function of the argument of perigee (Refer- 
ence 5-7). The effect of the long term trend on the perigee height of IMP 3 is 
shown in Figure 2. This long period trend rises for 1-1/2 years, reaching a max- 
imum of 41,000 km from the center of the earth; then it reverses and declines 
until perigee goes below the ground after a total lifetime of a little over 3 years. 

When a small portion of the lifetime is examined closely the changes in 
perigee height seem very irregular (Figure 3). During the lifetime, the dif- 
ference between one given perigee height and the following one may be as little 
as 20. km, or as much as 2200. km. These seemingly random fluctuations 
around the long period trend are due solely to the moon's gravitational field. 

They occur because the semimajor axis of the satellite is not negligible in com- 
parison with the moon's semimajor axis (a/a' - 1/3). The satellite's apogee 
moves only slowly in inertial space during the lifetime of the satellite, while the 
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moon is orbiting about the earth once every 28 days. Therefore when the satel- 
lite reaches apogee it may find the moon nearby; then on the next apogee pass 
6 days later, the moon will have moved nearly 90° away. It is this variation in 
the apogee-lunar configuration which is responsible for the erratic quality of 
the changes of perigee height, and, also, of all of the orbital elements. 

The variability of the lunar perturbation is of great consequence to the life- 
time near the beginning or near the end of the lifetime. Between launch and the 
first return to perigee, the perigee of IMP 3 is raised from 200 km to 600 km, 
bringing the satellite above the regions where the earth's atmosphere will have 
any effect on the orbit. 

Similarly, the erratic component in the lunar perturbation will determine 
on which orbit the satellite will be destroyed. At the 193rd perigee passage, the 
perigee height is 600 km. A very strong drop of 1100 km in perigee height oc- 
curs during the subsequent orbit, resulting in a predicted perigee height of 
500 km below the earth's surface, or, actually in the impact of the satellite be- 
fore perigee passage. 

The orbital period of IMP 3 also shows an erratic variation between 139 
hours and 143 hours (Table 1). This is ail unusual occurrence in eelcstiai 
mechanics; it happens because the satellite's trajectory extends out to 2/3 of 
the lunar distance. This phenomenon is predicted by numerical analysis, and 
is verified by the data. Because of the large and erratic variations, the term 
"orbital period" must be carefully defined; here, it is used to denote the time 
elapsed between one perigee passage and the next. "Perigee passage" is said 
to occur when the minimum value of the radius vector from the center of the 
earth is obtained. 

There is also a short period perturbation by the earth’s oblateness and 
higher harmonics. These perturbations are at a maximum at perigee and 
diminish as r -3 . Figure 4 shows a typical example of the short period oblate- 
ness perturbation on a prelaunch IMP orbit (Reference 2). The osculating ele- 
ments at perigee may predict a markedly larger period than occurs according 
to the definition of the period above. All of the orbital elements will show very 
sharp changes during perigee passage, however, most elements do not show a 
secular change following perigee passage. As the perigee is raised away from 
the earth, this effect diminishes. 

For most applications, therefore, it will be useful to find the minimum 
values of the radius vector. The difference between successive times of passage 
through the perigee should be used to compute the orbital period, rather than 
osculating elements, recognizing that this value will fluctuate 3-4 hours from 
one orbit to the next. 
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The inclination in a highly eccentric orbit can be perturbed through tens 
of degrees. The inclination of IMP 3 shows a long period trend, rising from 33° 
to over 50° (measured with respect to the plane of the earth’s equator) and de- 
creasing to less than 30° just before termination of the orbit (Figure 5). The 
tendency to rise rapidly initially and then to form a plateau is typical of highly 
eccentric orbits; the plateau may occur at very high inclinations. 

The latitude of perigee moves from 22°N to -10°S during the satellite’s 
lifetime (Figure (j). The longitude of the subsatellite point at perigee is random, 
because the satellite's orbital period is not commensurate with the rotation of 
the earth; the right ascension of the ascending node, which is fixed in inertial 
space, evolves gradually from +222° to +17(5° during the lifetime. 

It should be commented that there is no evidence of a secular change in the 
average value of the semi major axis. This remains an open question in celestial 
mechanics (Reference 8). Here the small apparent long term variation of the 
semimajor axis may l>e attributed to the earth's second harmonic coupled with 
the long-term evolution in perigee. 


IV. COMPUTATION 

A highly eccentric orbit provides a stringent test of a numerical integration 
program. Most major forces are required to compute the orbit accurately, along 
with careful computational techniques. 

The major gravitational forces are the sun, the moon and the earth, includ- 
ing oblatcness and higher order harmonics. The sun and the moon both perturb 
the perigee height strongly, and if either is inadvertently omitted, or incorrectly 
included, in a program, the perigee height will not be correctly computed. 

The oblateness of the earth has a small effect on the long term evolution 
of the orbital parameters. However, there is a very sharp short period effect 
of the oblateness perturbation which occurs very near perigee. Extreme care 
must be taken to prevent the development of computational inaccuracies during 
this portion of the orbit. Because the effect is so localized, the computation is 
very sensitive to the size of the integration interval near perigee. It can hap- 
pen that an interval which is satisfactory elsewhere in the orbit is far too large 
near perigee. In fact, most of the machine time necessary to compute an orbit 
is consumed during perigee passage. An improved method of calculating the 
oblateness perturbation at perigee would be very helpful in reducing machine 
time. 



If the integration interval is not controlled carefully, the orbit will receive 
an apparent (and erroneous!) impulsive force which will alter the orbit. After 
a few such perigee passes, the numerical accuracy of the computation will have 
become irretrievably degraded. The osculating value of the semi major axis is 
helpful in indicating the accuracy of the trajectory. In a correct computation, 
this value will rise sharply as perigee passage occurs, shortly after perigee, 
this quantity will drop back to the value shown just before perigee passage. 

This ’’spike” effect is to be expected, with the magnitude ot the spike depending 
on the height and position of perigee. 

If the integration interval is large compared with the spike, so that the 
spike is not computed correctly, the value of the semimajor axis after perigee 
passage will not return to the value before perigee. The effect is to add an in- 
correct impulsive change in the orbit. 

Large changes in the semimajor axis between two perigee passes are not 
necessarily due to numerical error. As was explained in the previous section 
on perturbation, the moon may impose a substantial variation on the semi major 
axis (that is, the orbital period). These variations should average out to near 
zero over a number of orbits. However, a drastic secular change in the semi- 
major axis should not occur if there is no substantial atmosphere drag and no 
passage around the moon occurs. For example, a decay of several earth radii 
in the semimajor axis would most likely indicate computational error. 

The computation of the ephemeris has been performed with the ITEM (Inter- 
planetary Trajectory Encke Method) program (Reference 10), using a Fortran IV - 
360 version. This program is based on a modified Encke method, which numer- 
ically integrates the perturbations on a reference ellipse. The reference body is 
the earth, with permissible gravitational perturbations including the moon, sun, 
Jupiter, Venus, Mars, the earth’s harmonics and the non- gravitational perturba- 
tion due to solar radiation pressure. 

There is an option in the program to employ a regularized variable as an 
independent variable rather than the time. The regularized variable is an 
advantage in computing highly eccentric orbits, since the velocity of the satel- 
lite varies by a factor of 40 between apogee and perigee. Integration intervals 
of equal time w ould allow' the satellite to go an excessive distance at perigee 
while barely changing at apogee. Therefore, a transformation to use as the 
independent variable is made: 


dt 

d/3 
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where t - time in seconds, f. the radius vector in km, and p = 398603. km 3 /sec 2 
(gravitational constant for the earth). Then, at perigee, taking p = 7000 km, a 
value of A/3 - 1/16 will equal a At = .6 sec; while at apogee,/? = 240,000 km, 
then for A/3 1/16, At - 25. sec. In other words, the use of A/3 will auto- 

matically contract the integration interval at perigee and expand it at apogee. 

The choice of A/5 - 1/16 yields the accuracy desired for this problem. This 
choice is indicated by the interval required to satisfactorily compute the earth's 
oblateness near perigee. A larger value of At could likely be used at apogee; 
but since most of the machine time is consumed performing the perigee computa- 
tion the value of At at apogee is not a significant factor. 

The value of A/j - 1/16 was compared with results from an otherwise identi- 
cal run utilizing A/3 - 1/32. The discrepancy in time of perigee passage between 
the two runs was an order of magnitude less than the differences between the 
computation and the data. Therefore, the value of A/3 = 1/16 was considered 
satisfactory. 


V. DATA AND TECHNIQUE 

Orbit improvement, rather than orbit determination, is used to obtain initial 
conditions which will produce an ephemeris valid over the entire two-year arc. 

A set of preliminary initial conditions are used to generate an ephemeris which 
is compared to the data and, from the differences between the data and the com- 
putation an improved set of intiial conditions is computed. 

The "data" itself is not ; ur observations, but rather is extracted from com- 
putations which are the result of previous reductions. That is, a sequence of 
short arc orbit determinations was available and the objective was to try to 
obtain a good fit to all the short arcs simultaneously, rather than to attempt to 
match the original raw data as obtained from the range and range-rate tracking. 
The short arcs are in the form of World Maps (Reference 9), which express the 
longitude ard latitude of the subsatellite point and the height above the earth's 
surface at intervals of one minute. These short arcs are determined over a 
period of about one month, that is, 5 or 6 satellite revolutions. 

Two quantities were selected from the World Maps for the purpose of com- 
paring the computation with the data: the time of perigee passage and the perigee 
height. These two parameters were of crucial significance in predicting the 
reentry of the orbit. Also, they provide a stringent test of the accuracy of the 
ephemeris. 

The time of perigee passage yields a critical evaluation of the ability of 
the ephemeris to compute an instantaneous position in space. If there is a 
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timing error in the ephemeris, for example originating from an incorrect esti- 
mation of the initial semi major axis, the computed value of the position in orbit 
will increasingly deviate from the actual value, although all other orbital ele- 
ments may be in good agreement with the "real” elements. 

The time of perigee passage provides a particularly useful parameter for 
evaluating the timing of the ephemeris because it is very well-defined in a 
highly eccentric ellipse and occurs when the best geometrical definition is ob- 
tained because the satellite is closest to the earth. In fact, it may be that the 
original orbit determination process will force the solution to be correct at 
this position by making compensatory errors elsewhere. Certainly, the ephemeris 
in this paper was obtained by a deliberate forcing of the solution to be accurate 
at the time of perigee passage. 

Since the World Map quantities are listed at intervals of one minute, an 
interpolation was found desirable to provide more accurate values of the perigee 
time and position. This allows a better and more convenient comparison with 
the numerical integration, which computes values at the instant of perigee 
passage. 

The height of perigee is also a well-defined parameter. It provides a meas- 
ure of the eccentricity and the variations of the eccentricity. It is critical to 
the lifetime of the satellite. Also, a good agreement in the height of perigee 
guarantees that the perturbing forces of the sun and the moon are correctly 
computed. The values of perigee were also obtained by interpolation from the 
World Maps. 

In addition, the latitude and longitude of the perigee position are compared. 
The latitude indicates the accuracy of the angular orbital elements. The longi- 
tude is a redundant measure of the timing error and the angular elements. The 
major usefulness of the longitude for a comparison program is to ensure that 
the rotation of the earth is included correctly. The initial orbital elements, 
obtained from the injection conditions as determined from the early orbits fol- 
lowing launch, were used to initiate a numerical integration of the orbits. The 
values of the time of perigee passage for every orbit were obtained from the 
computation and subtracted from the "observed" values. When these differences, 
or residuals, of the first 5 or 10 orbits were examined, a linear growth in the 
time of perigee passage would appear due to the accumulation of error when 
the orbital period assumed in the computation differs slightly from the actual 
period. The residuals between the observed and computed time of perigee pas- 
sage on the ith orbit after launch 

At = t - t 

p i p i obs p i comp 
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were summed by least squares to yield a change Ap in the orbital period P 0 
initially assumed 


Ap 


2t_ At 

p i comp p i 

it 2 

P • 

i comp 


p 


0 


which was converted into an improved value of the semimajor axis, a 0 ', by use 
of Kepler' s law 


# 



a 0 + Aa 


a o + 


2/3 Ap p- 


(time is measured from launch). The numerical computation was performed 
again, utilizing the new initial value of the semimajor axis, a 0 ' . This process 
is repeated until the residuals no longer exhibit a linear trend; then the com- 
putation may be extended over more orbital passages, and, if necessary, the 
linear trend may again be removed. Only a few iterations are necessary. It 
was found that by using 25 orbits, all linear accumulations were removed and 
that iterations through additional data would yield no further improvement. 

The values of the perigee height residuals were used to adjust the initial 
eccentricity. The major effect of changing the initial eccentricity is to move 
the base line of the residuals in perigee height up or down. 

It was found that a change in the initial value of the eccentricity, or of one 
of the angular elements (inclination, argument of perigee, or the right ascension 
of the ascending node) necessitated a change in the initial value of the semimajor 
axis although the value of the semimajor axis prior to the change in another 
orbital element had yielded satisfactory residuals in the time of perigee passage. 
This coupling of the semimajor axis with the other orbital elements is apparently 
due to the oblateness of the earth which strongly influences the energy of the 
orbit. (The oblateness causes a difference between the anomalistic and the 
osculating value of the period at perigee of about 4 hours on the IMP 3 orbital 
period of 140 hours.) However, satisfactory values of the semimajor axis can 
be obtained equally well after small changes in the orbital elements, providing 
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that first the remaining elements are modified, and then the iterations are per- 
formed with only an input change to the semimajor axis until the procedure 
converges. 

Due to this coupling of the corrections to the semimajor axis with correc- 
tions to other orbital elements, the initial conditions which provide satisfactory 
residuals in the time of perigee passage are not unique. For example, little 
change in the residuals is produced, if either set of orbital elements 


Set 1 Set 2 

a 0 21.726218 earth radii 21.726138 earth radii 

135273696 135272696 

e Q .95251083 Same 

i 0 33.828445 Same 

-138.45267 Same 

is used to initiate the numerical computation. It might be argued that this occurs 
because the time of perigee passage is only one observational parameter for 
testing the orbit and that consideration of two or more variables might provide 
a better fit of the orbit; however, this ’'observed” quantity is particularly sensi- 
tive and will in fact partially determine the choice of other quantities once a 
semimajor axis has been chosen. Further, not only corrections to the orbital 
elements but also model errors such as the value of the earth’s harmonics, 
may be eliminated by an appropriate choice of the semimajor axis. It should 
be stressed that the semimajor axis in Set 2 could not be substituted into Set 1 
and produce satisfactory residuals. A change of .5 km in a Set 1 would produce 
a substantial linear term in the residuals. The correlation between the semi- 
major axis and numerous other quantities constitutes a major difficulty in long- 
arc orbit analysis. 


VI. COMPARISON BETWEEN DATA AND COMPUTATION 

The residuals in the time of perigee passage provide the most critical evalu- 
ation of the accuracy of the ephemeris. The data from the first 25 orbits were 
used to adjust the initial value of the semimajor axis; then the numerical integra- 
tion was extended without further correction through for the 120 orbits during 
which there was radio coverage. The difference between the observed time of 
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perigee passage and the computed time for each orbit is shown in Figure 7. 

These residuals show no linear trend during the 120 orbits, indicating that 25 
orbits are sufficient to define the semimajor axis (or, more rigorously ex- 
pressed, to define the total energy). 

The residuals in the second portion of the 120 orbits exhibit a systematic 
deviation which resembles a low frequency sinusoidal. This deviation is cur- 
rently unexplained but may arise from one of the following sources: 

1. a small error in one of the initial conditions 

2. an error or omission in the model in the numerical integration, for ex- 
ample, an omission of a higher order harmonic in the earth’s gravita- 
tional field, or a poor choice of the value used in computing the solar 
radiation pressure force 

3. an error in the numerical integration. 

Of these choices, 2, is most likely. It can be difficult to locate the model error 
when the model is as complex as is required for this ephemeris. 

The largest residuals are no greater than 1 minute in time of perigee pas- 
sage, compared with an orbital period of 140 hours. This is an indication of the 
accuracy of the ephemeris. The systematic trend in the residuals shows the 
data itself is highly accurate. Since there is no discernible scatter in the data, 
the random error in the data cannot be evaluated, but it is certainly much less 
than one minute of time. It seems probable that this type of data is independent 
of systematic errors and that the systematic trend exhibited in the residuals is 
due to computation error. 

The residuals in the height of perigee are shown in Figure 8. Some estimate 
of how much of this deviation is due to the data can be made by examining over- 
lapping data arcs of 4 or 5 orbits; the residuals often change sharply when the 
arc is changed. Occasionally, there is a perigee passage which is computed 
during the overlapping position of two adjacent arcs. There are differences be- 
tween the overlapping arcs in the perigee height amounting to 6 km early in the 
lifetime and as large as 20 km toward the end of the data (when the perigee is 
near maximum distances from the earth); therefore, the data is at least that 
inaccurate. There appears to be a systematic trend as well; this is most likely 
due to a model error in the computation. The latitude of perigee residuals are 
less than .07°. 

Although the telemetry ceased operating after two years, there are visual 
observations confirming the latter part of the ephemeris. As the perigee height 
decreased, it became possible to acquire the satellite with Baker-Nunn cameras. 
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After almost one year without any type of observations, this ephemeris was used 
to predict the location of the satellite during perigee passage. One several 
occasions, the SAO tracking network acquired a satellite in the predicted loca- 
tion with the expected visual magnitude and with a spin rate compatible with the 
satellite’s characteristics (Figure 1). Although it is difficult to compare the 
observations precisely with the ephemeris, it is felt that the acquisition of the 
satellite implies that the time of perigee passage is correct to within ±3 minutes. 

The precision of the ephemeris and the fact that the last perigee before re- 
entry was 600 km above ground while the reentry "perigee” was 500 km below 
ground led to the hope that the reentry phase might be observed. A group on 
the island of Gan in the Indian Ocean attempted to observe the event (Reference 
11); however, there was a heavy cloud cover over the area that night and the re- 
entry was not observed. 


VD. EPHEMERIS 

The computed ephemeris is presented here for use as a standard in check- 
ing other programs. The ephemeris (Table II) is presented in two forms: 1) 
the osculating orbital elements are listed at perigee passage (to within 10" 4 
hours) for every 10th perigee passage during the satellite's lifetime: 2) the 
osculating orbital elements are listed every 30 days during the first year. (The 
tabulation is actually within a fraction of an hour of 30 clays. Usually the osculating 
elements change little during tMs time span.) The orbit number is defined by 
taking the "0th" orbit to be at launch. 

The geodetic longitude and latitude of the subsatellite point at perigee and 
the perigee distance in km from the center for the earth for the first 25 orbits 
are given in Table III. 

An ephemeris listing osculating elements at every perigee passage may be 
obtained from the author. 

A researcher interested in comparing his program with this ephemeris would 
most likely have one of two objectives: 

1. The ability of this program to compute the evolution of perigee height. 

If this is the researcher's objective, precision computation of the time 
of perigee passage is not necessary. The evolution of perigee can be 
computed by one of three methods: 

a. a special perturbations method such as employed here, which relies 
completely on numerical integration; 


12 


b. a ’’mixed" method which utilizes both analytic expressions and numeri- 
cal integration. An example of this is the Halphen method as employed 
by Musen (Reference 12), and applied by Smith (Reference 13). (This 
method computes only the long term evolution. Terms depending on 
the moon's longitude are omitted.) 

c. a strictly analytic method, such as developed by Lidov (Reference 6) 
as employed by Renard (References 1, 3). It should be pointed out 
that the IMP 3 orbit is a more stringent test than necessary for high 
eccentricity orbits with smaller semimajor axes. 

2. The ability of this program to track the satellite around its orbit. The 
accurate computation of time of perigee passage is a necessary, although 
not sufficient, condition to yield the correct position of the satellite. 

If there are model errors which are compensated for by an incorrect 
value of the semi major axis, there will be systematic errors in the 
residuals near, but not at, the perigee passage which are not necessarily 
evident in the computation of t . 

As has been discussed above, the data indicates that there probably are 
model errors in the generation of this ephemeris. The researcher may compen- 
sate for these and for model differences between his program and ITEM by 
adjusting the initial value of the semi major axis according to the procedure de- 
scribed above. 


VIII. CONCLUSIONS 

An ephemeris has been obtained which represents the "long" (monthly) arcs 
of a high eccentricity satellite with good accuracy. Probably with further 
study, and with the improved data now available on current satellites, even 
better precision can be obtained. 

The ephemeris can be used as a standard to evaluate a program's ability 
to compute either the long-term evolution of a high eccentricity orbit in the 
Earth- Moon system, or to track a satellite along its orbit. 

It may be possible to predict the longitude, latitude and time of reentry of 
future high eccentricity orbits far in advance of the event. The ability to predict 
will depend on the values of the final perigee heights. The IMP 3 was particularly 
favorable for prediction because of the large decrease in the perigee height dur- 
ing the final orbit. 
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Long arc determinations may be useful in evaluating the methods and the 
results of short arc determinations. In particular, it may be that ambiguities 
in the energy determination similarly affect short arc determinations so that 
equally good fitting initial conditions could be obtained with low errors ascribed 
to each but with larger differences between the two sets than the errors. Also, 
model errors may be compensated for by changes in the initial conditions pro- 
ducing an illusion of a good fit. It may be worthwhile to determine model errors 
by identifying observable quantities which are particularly sensitive to a model 
quantity— such as the application of resonant orbits to the earth’s harmonics. 

The time of perigee passage is an exceptionally sensitive parameter when 
it is available for large number of orbits and is therefore to be recommended 
for utilization in long arc studies. It insures that the instantaneous position of 
the satellite corresponds closely to the predicted position. Further, it appears 
that this quantity is determined very accurately when raw data has been reduced 
to produce a short arc of 5 or 6 orbits. The accuracy of the data is inferred be- 
cause there is no strong discontinuity in the residuals from short arc between 
one short arc and the next. 
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Table 1 


Variation of Time Between Perigee Passages 
Due to Random Lunar Perturbation 


Perigee Passage 

Time in Hours Since Last Passage 

1 

139.866 

2 

140.590 

3 

138.860 

4 

139.385 

5 

139.604 

6 

138.773 

7 

140.609 

8 

139.770 

9 

140.350 

10 

140.309 
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Table HI 


Perigee Condition* of IMP 3 


rbit f 

Longitude 

Latitude 

0 

-62.923 

21.564 

1 

-£*468 

42*923 

2 

38 • 60 C 

22.866 

3 

103*002 

42*2126 

4 

I7C.488 

22.772 

5 

-429*317 

22*632 

6 

-56.606 

22. 702 

7 

- 11*628 

22*342 

8 

46.018 

22.377 

a 

94*992 

-22*2 9 6. 

10 

144.894 

22.1 84 

n 

- 154*982 

22*354 

12 

-te.oi 6 

22. 1 68 

13 

-46*654 

22*156 

14 

-C.66S 

22.0C2. 

15 

46*720 

11*92 Pi 

16 

92.709 

21. 899 

17 

164*642 

21*656 

18 

- 132.246 

21. 752 

19 

- —72*906 

21.322 

20 

- 3.908 

21.51 1 

?1 

34.274 

20.957: 

22 

91.81 8 

2C.928 

23 

mo. te? 

2C.707i 

24 

- 170.118 

2C.477 

05 

-110.129 

20. 63 1 • 

26 

-62.626 

20. 1 46 

* T 






28 

21.661 

19. 85 C 

29 • 

47.623 

19,760 

30 

108.918 

19.636 

31 

174,324 

19^624 

32 

-12C.C11 

19.574 

34 

-4 1.749 


34 

7.7C8 

15.508 

35 

44r44f 

16.951 

36 

108.657 

19. 088 

37 

164^422- 

19,01 6 

38 

-140.016 

16.91 e 

39- 

-73.790 

to. i a o 

40 

-21.777 

16. 839 

44 

26.635 

16.91 5 

42 

72.008 

ie.764 

43 

114,035 

16,706 

44 

164. e«t 

16.59C 

48 

134.336 

18.574 

46 

-60.096 

18.480 

47 

—1.37* 

16.205 

48 

68.037 

ie.i69 

44 

97,431 

17,231 

80 

166.366 

17.204 


Perigee Distance 

6580.71 
6 05 a* 01 
7639.24 
7934*78 
791 6. 74 
0207*35 
8268.73 
9620 *<28 
971 7.20 

0441 -1? 

9903.49 

9704*73 

1 0545.22 
-1 0333. 99 

1 0406. 69 
1 Oft 0.3*66 

11274.44 
-l-151B*-69 

11510.26 
11621*87 

11994.70 
14236.42 

14374.23 

14616.92 
15117.77 

15266.92 
16903.51 

• I 7 0 4 7 . 03 
1 7467. 76 

1 7981 .45 
1 861 0.51 
1 9 22 4.61 
1 934 0. 30 

20005.67 

22246.86 

22306.37 

22345.00 

22686.98 

32656.60 
24077.20 
23867.32 
23926.82 

24147.27 

24248.45 
246 7 0.02 
24575.04 

2482 1.60 
25001.31 
272 32. 71 
27358.16 
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Figure 1. Photograph of IMP-III taken by Smithsonian Astrophysical Observing Station at Olifantsfontein 

South Africa, on March 27, 1968 at 18:24:26-5 UT 
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Figure 2. Evolution of Perigee Height During 3-Year Lifetime 
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Figure 3. Random Effect of Short-Period Perturbations on Perigee Height 
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Figure 4. Perigee Perturbation of the Semimajor 
Axis of on IMP-Type Satellite 
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INCLINATION IN DEGREES 



Figure 5. Evolution of Inclination of IMP III 
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LATITUDE OF PERIGEE IN DEGREES 



Figure 6< Latitude of Perigee During Lifetime of IMP 3. 
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Figure 8> Reeidual* in Perigee Height 
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